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The flavor-changing electromagnetic dipole operator O7 gives the dominant contribution to the 
B — > X a j decay rate. We calculate two-loop QCD corrections to its matrix element together with 
the corresponding bremsstrahlung contributions. The optical theorem is applied, and the relevant 
imaginary parts of three-loop diagrams are computed following the lines of our recent t — > XbW 
calculation. The complete result allows us to test the validity of the naive non-abelianization 
(NNA) approximation that has been previously applied to estimate the NNLO QCD correction 
to F(B — » X s j)/F(B — + X u ev). When both decay widths are normalized to ml R in the same 

■ renormalization scheme R, the calculated 0(o? s ) correction is sizeable (~ 6%), and the NNA estimate 
\^ ' is about 1/3 too large. On the other hand, when the ratio of the decay widths is written as 

S x m 2 b —(mb)/m1 pole , the calculated 0{o? a ) correction to S is at the level of 1% for both the 

■ complete and the NNA results. 
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PACS numbers: 13.40.Hq,14.65.Fy,12.38.Bx 



I. INTRODUCTION 



The radiative decay b — > occurs only through quantum loop effects, similarly to the muon anomalous magnetic 
moment (g M — 2) or radiative hyperon decays 0. The largest contribution by far to (g^ — 2) is due to electromagnetic 
interactions. In contrast, the flavor-changing hyperon (s — > d) and B meson (b — > s) decay amplitudes involve heavy 
particles such as the W boson and are consequently very rare. In particular, the b — > 57 transition can be predicted in 
the Standard Model with good accuracy, and it offers a relatively low-background probe of possible new phenomena 
such as supersymmetry (e.g., see Ref. 0|)- 

High-accuracy measurements of the B — > X s ~/ rate in B-factories Q warrant sophisticated calculations of high-order 
Standard Model contributions. Electroweak loop effects, without which this decay would not occur, have now been 
studied to two-loop accuracy 0, H, IE EI ■ F° r the precise determination of the rate, the QCD effects are crucial and 
are fully known to the next-to-leading order (e.g., see Refs. I9|. [loll. Potentially important effects resulting from the 
binding of the b quark in the B meson were explored in Ref. [Ill ] . 

At present, several groups are working at the determination of the next-to-next-to-leading order (NNLO) QCD 
corrections. A discussion of the various required studies can be found in Ref. Since that review was published, 
new ingredients have been provided: all the relevant three- loop anomalous dimensions \lA Il4j. three-loop matching 
conditions [l5l | , as well as certain counterterm contributions to the three-loop matrix elements of four-quark operators 

E3 ■ 

Among the most challenging missing quantities are the two- and three-loop matrix elements of several operators. So 
far, the only complete two- loop O (cv s ) results exist for the four-quark operators, where one of the loops is a fermion 
loop H^l Ell- I n addition, parts of two- and three-loop O (oif) matrix elements, involving gluon vacuum polarization, 
were found in Ref. |l9j | . This last result is very important since it automatically determines corrections to the matrix 
elements of order ot 2 s (3o, where /3o = 11 — 2Nf/3 is a large parameter (Nf = 5 denotes the number of quark species). 

In the present paper, we provide the first calculation of full two-loop corrections to the matrix element of the 
electromagnetic dipole operator O7 that is responsible for the lowest-order b — > sj decay rate. We reproduce the 
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FIG. 1: Three loop self-energy diagrams required for the O (q?) computation. Thick solid closed loops in (ab,ac,ad) depict 
massive (b quark) loops. Thin solid closed loops in (y), (z), (aa) denote massless fermions. From the latter diagrams we find 
also the contributions of gluon and ghost loops, as explained in the text. 

a 2 s fio effect found in Ref. [l^ and also provide the a 2 s corrections that are not enhanced by /3q. With this result, we 
can check the extent to which the (3q effect is dominant. The conclusion turns out to depend very much on what 
renormalization scheme is used for the overall factor of mjj in the expression for the decay rate. 

Since the result obtained here is only a partial contribution to the future full NNLO correction, we present a 
number of intermediate results and describe in detail our renormalization procedure. We hope that this will simplify 
the utilization of our result when the remaining ingredients are known at the NNLO level. 

The paper is organized as follows. In Section^] we describe the calculation of the relevant diagrams as well as their 
renormalization. Section [IIII is devoted to presenting our main results and discussing the large-/?o approximation. We 
conclude in Section llVI Contributions from particular diagrams are listed in the Appendix. 
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II. THE CALCULATION 



The decay rate of b into s, 7, and up to two gluons or a light quark pair can be written as 



r(6 - xr ton 7K>£o 



G 



p u em 



32tt 4 



KI 2 E^ (/*) Cf 00 G y -(i?o,M) 



(1) 



where mj, is the pole-mass and mb{u,) is the MS running mass of the 6 quark. The effective [2G| Wilson coefficients 
in the relevant low-energy theory are denoted by C? ff (/i). The photon energy cutoff E is assumed to be significantly 
below the endpoint, i.e. mj, — 2Eq A QCD . This is a necessary condition for the perturbative decay width Q to be 
a good approximation for T(B — > X s j)e^>e - 

In this paper, we focus on the contribution of the operator 



Or 



to the decay rate. More specifically, we calculate 



G 77 (0, M ) = 1 



4tt 



A, + ( - ) A 2 
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Our result for (777(0, u) can be combined with the very recent findings of Ref. [2l| to obtain G^{Eq, u) at the NNLO 
for any value of Eq that is sufficiently far from the endpoint. 

As far as the coefficient X\ is concerned, we confirm the well-known result of Ref. j^. The NNLO correction X 2 
can be subdivided into color structures, 



X 2 = Cf (T r N l X l + T R N H X H + CfXa + CaXna) 



(4) 



where Cf = 4/3, Ca = 3, and Tr = 1/2 are the SU(3) color factors, while Nl and Njj denote the number of light 
(m q = 0) and heavy (m q = nib) quark species (Nh + Nl = Nf). 

Let us now briefly outline the method applied to the calculation of X 2 . We use the optical theorem to map all the 
virtual corrections and real radiation contributions onto a system of self-energy diagrams as described in Ref. \'2 i in 
the context of the decays t — » XbW and b — ► X u lD. The topologies which have to be taken into account turn out to 
be identical to those considered there for the top quark decay. 

We consider only gluons of virtualities of order mj. In the effective theory, this is the only mass scale. In other 
words, we do not consider hard gluons that would resolve the structure of the effective vertex sbj since their effects 
are accounted for in the Wilson coefficients C 4 off (/x). Diagrams needed for the present calculation are analogous to 
those of the hard asymptotic region for the top quark decay, studied recently in Ref. |U 0] ■ 

Diagrams contributing to X2 are presented in Fig. ^ All the particles except for the b quarks are treated as 
massless. Both the UV and the IR divergences are regulated dimensionally in D = 4 — 2e dimensions. The results 
for /i = nib are collected in Table I in the Appendix along with the relevant color factors. They have been computed 
in a general covariant gauge. The Feynman gauge results can be obtained by setting £ = 0. Diagrams that are not 
symmetric under left-right reflection are already multiplied by 2. The quantities ln4-7r and 7 that account for the 
difference between the MS and MS schemes are omitted in Table I and in the remainder of this section. 

Diagrams with closed gluon or ghost loops are not shown explicitly in Fig.^ since their contribution can be found 
from the total contribution of the light fermions. To this end, after computing the light-fermion diagrams, we replace 
(see for example Ref. 0) 
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Let -B 3 denote the sum of all the three-loop diagrams from Table I, after performing the above replacement. Our 
final result for 6*77(0, u) is found according to the following formula 



where 
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stand for the one-loop (Born) diagram and the sum of two-loop diagrams, respectively. The MS rcnormalization 
constant for the operator vertex 



Z MS Z MS _ 1 , C F /Q s \ 2 Cf 

m 77 4jt £ \4irJ e 



257 19 13 1 /l 11 2 
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is found from the anomalous dimensions published in Ref. 26]. The on-shell rcnormalization constant of the quark 
field can be written as 1271 
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where 
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(13) 
(14) 
(15) 

+ 0(e). (16) 



Mass rcnormalization in the 6 quark propagators is accounted for by squaring these propagators in the two-loop 
diagrams, which turns B 2 into 



B™ = C F 
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(17) 



In the expression JfjJ) for G77, the above quantity gets multiplied by Z,° s — 1 = j^Pi + O (a: 2 ,). For completeness, the 
one-loop gauge coupling renormalization constant should also be mentioned 
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III. RESULTS 



Our final results for the contributions to 6*77(0, fx) read 

X X = Cf( 1 4+^-^ 
\ 3 \i 3 
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The complete (logarithmic and constant) contribution of the light quark loops has already been found in Ref. |19j . 
However, the decay width was normalized there with mjj rather than with m b (fi) 2 ml as in Eq. Q here. In order to 



compare with that study, we multiply our result for G77 by m 2 (/z)/m 2 . In other words, we write 

G 2 P a em ml 
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The connection between the pole-mass m;, and the MS mass m b (n) is now known to the three-loop order [28| . Here, 
we only need it to two-loops [2!j 
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Using the above relation, we obtain 
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We find complete agreement of the X^ result with Ref. [13]. In that work, along the hypothesis of naive non- 
abelianization (NNA), Xl was multiplied by —3/2 P$(Nj = 5) in order to estimate X 2 - Our complete result for X 2 
allows us to check this hypothesis. Including all the SU(3) color factors, our analytic result leads to (for Nl = 4 and 
N H = 1) 
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(25) 
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Evidently, there are substantial differences between these expressions, from which we conclude that the NNA hypoth- 
esis does not necessarily improve on the Xl component of the O (a 2 ) part of the calculation. 

The large numerical value of the NNLO correction coefficient in Eq. 1)25(1 may be traced back to the infrared 
sensitivity of the pole mass m b whose fifth power stands in front of the expression ((20(1 for the decay rate. Following 
Ref. [33, we shall normalize the b — ► XP arton 7 rate to the semileptonic rate 
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From the results Ref. |3lJ, one finds 

G u ~ 1-9.65 
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Dividing our results by T(b — ► XP arton ei/) and expanding up to (a 2 ), we obtain 
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Note that (G77/G m )nna differs from (G77)nna/(G m )nna- The latter quantity gives a worse approximation to the 
complete result. The same is true for G77. 

In order to indicate where the renormalization of mb matters in the above expressions, we have introduced the 
superscript "pole" for the pole mass. Actually, no renormalization of needs to be performed when evaluating the 
O (a 2 ) terms in the NNA approach. However, it is mandatory to identify the mass in this approach with the pole 
mass because it often originates from the square of the external momentum. 

Comparing the /it-independent terms in Eqs. (|31|I - H34|I one concludes that the perturbation series converges much 
better and the NNA gives a better approximation for GnjG u rather than for Gjj/G u . This observation confirms 
that the normalization of the top quark contribution to the b — > S7 amplitude which was applied at the NLO in 
Ref. indeed helped in reducing the NNLO contributions that were unknown at that time. As far as the charm- 
sector amplitude is concerned, no conclusion can be drawn yet, because several important NNLO ingredients are still 
missing. 

We have checked that the /z-dependent terms in the complete (i.e., non-NNA) expressions for 677(0, /it) = 
677(0, ^t)m 2 (/x)/m 2 cancel out (analytically) with the corresponding ones that originate from the Wilson coefficient 
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Of course, it does not mean that the quantity (Gf ff (^t)) 2 G77(0, fx) is ^-independent at O 
need to be included for a complete cancellation, for instance, 
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As far far as the /^-independent terms in the Wilson coefficients are concerned, wc can check their values for the 
top-sector amplitude, for which all the relevant Wilson coefficients are now available at the NNLO [13, E3, 0- In 



7 



particular, setting the matching scale /io to m t (m t ) and the low-energy scale fi to m b (m b ) , we find 
C^ ff (m b (m b )) = Cf )cS (m b (m b )) 

Cl cS {m b {m b )) = Cf )0ff (m fc (m fc )) 

Thus, no large corrections to the Wilson coefficients are being observed, which means that the NNLO QCD corrections 
to (C£ cff ) 2 G 77 /G u arc significantly smaller than to (C£ off ) 2 G 77 /G u . 

Among the dipole operator contributions, there are still missing two-loop matrix elements of Os and the O (a 2 ) 
corrections to the interference of amplitudes arising from 7 and Os- The set of master integrals in the form available 
so far [24| is not sufficient for those calculations. The reason is that the imaginary parts of those integrals are presented 
as sums over all cuts, while in the case of Os we sometimes have cuts which do not correspond to the decay b — > S7. 
Thus, it would be desirable to recalculate the master integrals in such a way that each individual cut contribution 
is known separately. If that were done, one could apply the same algebraic reduction of all integrals to the set of 
master integrals, keeping track of the relevant cuts. This would give an analytic result for G 7 s(0,[i). A calculation 
of Gs8(Eq,/j,) would be much more difficult because of the IR divergences at E$ — > and collincar divergences at 
m a — > 0. Fortunately, the effect of Gss on the decay rate is suppressed by the square of the down quark charge or, 
more precisely, by (QdCs/C 7 ) 2 ~ 0.03. Consequently, the O (a 2 ) corrections to Gss(E ,fi) are negligible. 

IV. CONCLUSIONS 

We have evaluated two-loop QCD corrections to the matrix element of 7 together with the corresponding 
bremsstrahlung contributions. The size of the resulting (partial) 0(a 2 s ) correction to T(B — > X s j)/T(B — > X u ev). 
depends very much on the conventions for the factors of m b that normalize the decay rates. When both of them are 
normalized to ml R in the same renormalization scheme R, the 0(a 2 s ) correction is sizeable (~ 6%), and the NNA esti- 
mate is about 1/3 too large. On the other hand, when the ratio of the decay widths is written asSxra^ —[m b ) / m\ pole , 
the calculated 0(a 2 s ) correction to S is at the level of 1% for both the complete and the NNA results. 
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APPENDIX A: RESULTS FOR PARTICULAR DIAGRAMS 



Diagram 


Color factor 




(a) 




\ [32 32? + 8? 2 ] + 1 [ 1264 10 fC+ 2 ^ 2 ] + T 32tt 2 ^ 167T 2 ? 2 40^ 2 14 f 4 ? + 41 9 56 ? 2 


(b) 




-92 + 192C 3 ? + 2^C3 - 48tt 2 ? + SMvr 2 In 2 - 8tt 2 - i^vr 4 + 96? + 16? 2 


(c) 


c F (c F - 


jf [16? + 16? 2 J + i [64 - ^vr 2 ? - ^tt 2 + if 6 -? + ^? 2 J + if 2 - - 448?3? - 448?3 - ^vr 2 ? 
- 32tt 2 ? 2 - 144tt 2 + if 7T 4 + Ifl? + if2? 2 


(d) 


c F (c F - at) 


4[ 16 32? 16? 2 ] + 4 [ 296 6 fe 3 fC 2 ] 47 9 24 +64^ 2 ? + 32^ 2 ? 2 +32^ 2 10636 ? 53 9 18 ? 2 


(e) 


cl 


$ [16? 2 ] + - [- f ^ 2 ? + 32? + 2§°? 2 ] +128-256C3?- ^ ^ 2 ?-32^ 2 ? 2 - f tt 2 + f|^ + 272?+ ^? 2 


(f) 


c F (c F - %) 


i [_ 64 + 64? - 16? 2 ] + 2« - 512(3 + 96^ 2 ? - ifn 2 - ^? - ^? 2 


(g) 


cl 


4 [16 + 32? + 16? 2 ] + 1 [ 320 + 640 ? + 320 ? 2 ] + 5024 64^ 2 ? 32^ 2 ? 2 32^ 2 + 10048 ? + 5024 ? 2 


00 


cl 


£ [8 + 16? + 8? 2 ] + \ [i|2 + 2M? + if 2 -? 2 ] + 2S5 _ 327r 2 £ - l 6 ^ 2 ? 2 - 16^ 2 + ^? + ^f 8 -? 2 


(i) 


cl 


\\ 32 16? + 16? 2 ] + 1 [ s f 536 ?+ 320 ? 2 ] 16 ° 12 +32tt 2 ? 327r 2 ? 2 +647r 2 10844 ?+ 51 g 68 ? 2 


0) 


cl 


4 [64? - 32? 2 J + i [64 + ^tt 2 ? - I^tt 2 + 1^6? - ^? 2 J + 832 + 640^3? - 1280(3 + ™ tt 2 ? 

+ 48^ 2 ? 2 896 tt 2 + 26840 ? 10192 ? 2 


(k) 


Cf (C F - at) 


-560 - 576Cs? + 2304? 3 + 16tt 2 ? 2 - ^f-Tv 2 + § tt 4 + 208? - 16? 2 


(I) 


Cl 

^ F 


[-16? - 16? 2 ] + \ [-32 + f tt 2 ? + f tt 2 - ii 6 -? - 2|2? 2 ] - 352 + 256?3? + 256?3 + ^tt 2 ? 
+ 327T 2 ? 2 + 22° tt 2 - Sf£? - Sip? 2 


(m) 


cl 


jz [-32? - 32? 2 J + \ [-32 + ^tt 2 ? + ^tt 2 - -^f 6 - ? - 2|2? 2 J - 272 + 256(3? + 256(3 + ^tt 2 ? 

+ 64tt 2 ? 2 + 320 tt 2 12496 ? 10048 ? 2 


w 


cl 


jj [16? 2 ] +i [-f ^ 2 ? + 32? + 2|2? 2 ] _l28-256C3?-^^ 2 ?-f ^ 2 ? 2 + if ^ 2 -l^ 4 +368?+2f°? 2 


(o) 


C F (C F - %t) 


i [64 + 16? 2 ] + ifS + 192?3? - ^Cs - 48^ 2 ? - f tt 2 ? 2 + if 8 -^ 2 In 2 - 224^ 2 + 2M^ 4 + if? 2 



TABLE I: Imaginary parts of three-loop self-energy diagrams in a general covariant gauge. 
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Topology 


Color factor 




(P) 




\ [736 400? + 16c; 2 ] + 1 [ 86 3 72 32 3 48 ? + 320 ? 2 ] + 136 9 640 + 400tt 2 ? 16tt 2 ? 2 736tt 2 62 f 8 ? 

, 5024 pi 
" r 9 ? 


(q) 




\ [128 32? 16? 2 ] + 1 [ 21 3 28 + 224 ? 320 ? 2 ] + 14264 2 48Ca + 32tt 2 ? + 16rr 2 ? 2 128tt 2 

1 12416 c 5024 c 2 
+ 9 ? 9 ? 


(r) 


cl 


[-576 + 272? - 16? 2 ] + \ [-1280 + f rr 2 ? - ^rr 2 + if 8 -? - 2§2? 2 ] - 7184 + 256Cs? 
2048C3 5 f ^ 2 ? + 16tt 2 ? 2 + 320 tv 2 + 35296 ? 5024 ? 2 


M 


C F (C F - 2a) 


4 [16? + 16? 2 ] + 1 [64 64 7T 2 ? 64 rr 2 + 446 ? + 320 ? 2 ] 20680 256? 3 ? + 384? 3 566 tt 2 ? 
- 16rr 2 ? 2 + 384rr 2 In 2 - 2§* rr 2 + f ^ 4 + ^f 2 ? + 5fi? 2 


(t) 


— \C f Ca 


^ [-48? - 24? 2 ] + i [-96 + 32tt 2 ? + 64tt 2 - 432? - 176? 2 ] - 1040 + 672?3? + 768Ca + 256vr 2 ? 
+ fn 2 e + ifV - f 2 -^ 4 - 2f°-? - 968? 2 


(u) 


— \G f Ca 


272 + 288Cs? - 192Cs + ^tt 2 ? - 16rr 2 ? 2 + 88rr 2 - ^tt 4 - 320? 


(v) 


— \C f Ca 


\ [ 192 + 48? + 24? 2 ] + 1 [ 1424 8? + 192? 2 ] 24266 52( V? 32tt 2 ? 2 + 848 rr 2 60?+ 35 3 84 ? 2 


(w) 


— \C f Ca 


[48 + 72? + 24? 2 ] + 1 [408 + 540? + 168? 2 ] + 2f°- - 96Cs? - 192^3 - 144rr 2 ? - 48tt 2 ? 2 - 96tt 2 

+ 3058? + 22p? 2 


(x) 


—\CfCa 


£ [-48? - 24? 2 ] + ± [-96 + 32rr 2 ? + 64tt 2 - 432? - 176? 2 ] - 960 + 480< 3 ? + 1152< 3 + ^tt 2 ? 

+ If.**? + 4327T 2 - 2f 7T 4 _ 8008^ _ ggg^J 


(y) 


TrC f 


£ [-32] + i [-352] - 25M + f 


(*) 


TrC f 


h m + ¥ 


(aa) 


TrC f 


7* [If] + 7 [ U 9 4 + "fA + 12 2 9 7 8 + 17 3 92 C3 + IfTT 2 


(Bb) 


TrC f 


^[_ 64 ] + I[-ifi] _21 f 8 +64?r 2 


(ac) 


TrCf 


i[16] + § 


(ad) 


TrC f 


^ [ f ] + I [ lOjii + 256 ,2] + 3^ + ^ _ 64^2 



TABLE I: Imaginary parts of three-loop self-energy diagrams in a general covariant gauge (cont). 



